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00 ' Abstract 

A particle that moves along a smooth track in a vertical plane is influenced 
C^ | by two forces: gravity and normal force. The force experienced by roller 

<**> I- coaster riders is the normal force, so a natural question to ask is: what shape 

of the track gives a normal force of constant magnitude? Here we solve this 
problem. It turns out that the solution is related to the Kepler problem; the 



trajectories in velocity space are conic sections. 



1 Introduction 

Motion along curved trajectories is always associated with a normal force which 
according to kinematics is given by, 

v 2 
F n =m—, (1) 

P 



where v is the speed and p the radius of curvature of the track. If the only physical 
force is the constraint (or reaction) force from the track, constant normal force 
simply requires v ex ^fp. 

In most cases, however, the problem of interest involves a conservative force, 
usually gravity, and a normal (constraint) force. In these cases the problem is 
conservative since the normal force does not do work, by definition. One applica- 
tion is to roller coasters and these have recently been discussed in the pedagogical 
literature, see Pendrill pQ and Miiller [2J. In simple mechanical problems given to 
students the shape of the track is often taken to be circular but in practice this leads 
to unpleasantly large time variation of the normal force. In practice therefore more 
complicated shapes are used, in particular the clothoid since the curvature (XI P) of 
this curve varies linearly with arc length [5J [3] ■ Students that wish to understand 
the physics behind amusement park experiences in greater detail must therefore be 
prepared to study curves more advanced than the circle. 

Here we will derive the shape of the tracks that produce a constant normal 
force. Pendrill [1 formulates the problem but does not solve it. Auelmann [4] 
has studied the effect of adding a normal force of constant magnitude to the central 
gravitational force of the Kepler problem. He showed that this problem is integrable 
and derived some general results. Adding a normal force in space requires the firing 
of rockets perpendicular to the velocity. The much more natural problem of a 
constant gravitational force, a good approximation in amusement parks, with the 
normal force coming from a smooth track, also turns out to be integrable, as we 
will show in this article. 
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2 The problem and its conserved quantities 

Gonzalez- Villaneuva et al. [5] found that the trajectories of the Kepler problem are 
arcs of circles in velocity space. As it turns out the trajectories of our problem are 
circles in acceleration space. Integrating once we thus find that the trajectories in 
velocity space correspond to the conic section solutions of the Kepler problem. The 
shape of the tracks that lead to a normal force of constant magnitude N are then 
obtained by one further integration. 
We write the equation of motion, 

mi) = —mge y + N eg, (2) 

where eg = — sin 9 e x + cos 9 e y , and 9 is the angle between the velocity vector and 
the x-axis. We use polar coordinates in velocity space and put v — x e x + y e y = 
\/x 2 + y 2 (cos9e x + sin9e y ) = re r . Clearly the vector eg is perpendicular to the 
velocity r e r and thus normal to the trajectory. 

We note that ([2]) means that the tip of mi) is on a circle of radius N centered 
on the tip of the gravitational force vector —mge y . The two Cartesian components 
of this vector equation are, 

rait — —Ny/\fx 2 + y 2 , my = —mg + Nx/ 'y x 2 + y 2 . (3) 

We introduce the following dimensionless variables, 

Q Q Q XV 

x <- -»j;, yi ~y, dt < dt, v x i , v y < , (4) 

Vq Vq v v Vq 

and put, 

A = N/mg. (5) 

Using this the equations of motion are, 

v x = -Xvy/^jvl+vl, v y = -l + \v x /^Jvl + v 2 y . (6) 

We denote the dimensionless speed by r so that, 

v x = rcos9, v y — rsin#. (7) 

Transforming Eqs. ^ to velocity space polar coordinates r, 9 gives, 

r = -sin0, r9 = X-cos9, (8) 

which means that, 

dr = -sin6»dt, r d.6 = (A - cos 9) dt (9) 

Add the first of these multiplied with A — cos 9 to the second multiplied with sin 9 
to get, 

(A - cos 6) dr + sin 9 r d9 = => d [r (A - cos 9)] = 0. (10) 

Using the second of the Eqs. © then shows that, 

d[r(A-cos0)] =d(r 2 9) = 0. (11) 

Evidently, 

r (A -cos (9) =r 2 9 = L, (12) 

is a conserved quantity (constant of the motion) . One recognizes this as a quantity 
that corresponds to conserved angular momentum (or sectorial velocity) in velocity 
space. 



Obviously there is another conserved quantity for this problem: the energy 

Without loss of generality we can assume initial conditions such that E — 0. We 
then have, 

r 2 
y=-y (14) 

We thus have a two dimensional problem with two conserved quantities and this 
means that we have an integrable problem. 



3 Velocity space trajectories 

From Eq. (dD we find that, 

L/X 



1 - (l/A)cos( 



(15) 



r and 9 being polar coordinates for the velocity vector. This is easily recognized 
as the expression for a conic section (see e.g. Goldstein [B]) where 1/A corresponds 
to the eccentricity, and where the origin corresponds to one of the foci. We have 
arrived at the conclusion that the velocity space trajectory, the hodograph, of a 
constant normal force trajectory is a conic section, i.e. an ellipse, a parabola, or a 
hyperbola. 

We note that, only in the elliptic case, when the eccentricity 1/A < 1 will the 
speed, r = s, be finite for all angles. For A = 1, the parabolic case, the speed goes to 
infinity for cos — > 1. When A < 1 the speed goes to infinity along the asymptotes 
9 — > ± arccos A. 



4 Integrating the problem 

Since this two dimensional problem has two constants of the motion it is fully 
integrable. Here we use this to calculate the trajectories for the various cases that 
arise. We chose initial conditions so that, at time t = 0, we have, ro = 1, and 
xq = 0. According to fHI) we then also have j/o = — r^/2 = —1/2. 

4.1 The trivial straight line case 

The simplest case is the trivial case of L = 0, or A — cos9q = 0. From (fT2|) we then 
have that 9 = so 9 is constant = 9q. Then f = — sin#o and we have, 

= O , r = l-sin9 t. (16) 

Using that r — s, where s is arc length, we find, 

s(t) = *-sin0 o -, x(t) = cose (t-ahx9 -j, y(t) = - ^— ^M . (17) 

The second of these follow from time integration of v x = x = tcos9q, and the 
second from (fTi)) . These results can also be written, 

x(s) — cos0o s, y(s) = — 1/2 + sin0 o s, y(x) — ta,n9 x — 1/2, (18) 



so we are dealing with a straight line of slope tan 6* = Vl — A 2 /A. 



4.2 The non-trivial cases 

If L = A - cos O ^0, then 

A-cos0 o 

r = ~ x H"' ( 19 

A - cos 

according to Eq. (fi~2")) . Using this and the second of Eqs. © one finds that rd0 = 
(A — cos9)dt and hence that, 

A -cos 0q (A-cos0 o ) 2 (A - cos0 o ) 2 cos0 

(A — cosfc*)^ (A — cosfl)-* (\ — cos9y 

The second of these is obtained by noting that ds/dt — r, and the third follows 
from dx/dt — r cos 9. 

Finding the functions t(9, A), s(0, A), x(8, A) thus requires the following integrals 
to be done: 

( = <a —«I(a^¥ < 21 > 

^(A-costf^ VtLy (22) 

„ x2 Z" 9 cos0'd0' 
- (A - COS ^X (A-cos^ - (23) 

A complete expression for the trajectory is obtained by noting that the function 
2/(0, A) is given by, 

y = ^- cme ^ 2 2(x-Ler (24) 

according to Eqs. ([13]) and (HT)|) . 

When doing the integrals (|2T|) - ([23)) it is convenient to distinguish the following 
cases: 



Case 1 
Case 2 
Case 3 

Case 4 
Case 5 



A > 1, -7T < < 7T, 

A > 1, < < 2tt, 
A = 1, < < 2tt, 
< A < 1, — arccos A < < arccos A, 
< A < 1, arccos A < 9 < 2n — arccos A, 



o=0 


(25) 


= 7T 


(26) 


= 7T 


(27) 


0o=O 


(28) 


O =7T 


(29) 



The explicit results are given in Appendix A. 

4.3 The elliptic case 

It turns out that the complicated expressions obtained by direct integration using 
the variable become much simpler if one introduces a suitably chosen new inte- 
gration variable. In the Kepler problem this new variable is called the eccentric 
anomaly [BJ[7]. 

For Case 1 we chose, 



/ /A + l sin0 \ . . 

^ = 2arctan(y x - TTT -^, (30) 

so that, — 7r < ip < it, is in the same range as 0. One then finds that e.g. dt, as 
given in Eqs. (|20p. becomes, 

dt=(A-cos0 o ) ^ 2 "_ C °^ d^. (31) 




Figure 1: Plot of Case 1 trajectories for A = N/mg = 2 (red curve with two loops) 
and for A = 4 (green curve with four loops). As in all our plots the a;-axis is 
horizontal and the y-axis vertically upwards. 



Also the other integrals simplify considerably. 

For Case 1 we then get the following explicit results for the trajectory as a 
function of the parameter tp, 



x(tp 7 X) 

yftM) 



(A-l) 



(Xip + sin^>) 



(32) 
(33) 



(A 2 - l) 3 / 2 

affll^ [( 2A ' + W+ (4A + cosV)sin^ , 

2(A^- l) 5 /2 i 3A ^ + l 2 ^ 2 + *) + Xcoa ^ sin ^} > ( 34 ) 

(A-l) 2 



2(A 2 - 1) 



(A + cosVO i 



(35) 



from Eqs. (|21|l - (|24[) . One notes that one can allow ip to take values from — oo to 
oo here. 

Two such curves are plotted in Fig. [T] One notes the limit A — > 1+ which cor- 
responds to a straight horizontal line. For this case therefore x(i/j, 1+) = s(ijj, 1+). 
In the limit A — > oo the curve approaches a circle. 

Let us get some explicit numbers out. Using the factors given in Eq. Q, vq 
being the speed at the bottom of the loop, we now find that the time T that it 
takes to traverse a complete period is given by, 

(36) 



Vq, 



t'o 



9 9 VX 2 - 1(A + 1) 

where the function t(ip, A) is given by Eq. (|52"j) above. Similarly the total length L 
of one period of the track is given by, 



9 9 VX 2 - 1(A + 1) 2 



the horizontal width W is, 



3Att 



W=^[x(2n,X)-x(0,X)] = ^-^-^ , 

9 9 VX 2 - 1(A + 1) 2 



and the vertical height H is, 



2A 



H=^,X)-y(0,X )] = f {x + ir 



(37) 



(38) 



(39) 



We also find from (fK)|) that v(8) = Vq{\ — cos6*o)/(A — cos#), so the speed v v at the 
top of the loop [8 — it) is given by, 

(A-l) fAn . 

v ^ = Va lxTTy (40) 

Let us take the speed at the bottom of the track to be Vo — 20 m/s (this corresponds 
to 72 km/h) and use g = 9.81 m/s 2 . One then finds that, for A = 2, i.e. normal force 
N = 2mg, T = 4.93 s, L = 73.96 m, W = 49.30 m, # = 18.12 m, and -^ = 6.67 m/s. 
For A = 4 this becomes T = 2.65 s, L = 43.66 m, W = 15.88 m, H = 13.05 m, and 
v w = 12.00 m/s. 

Finally we note that Case 2, the other elliptic case with A > 1, simply gives 
similar displaced trajectories compared to Case 1. For Case 2 the top of a loop will 
be at, x = 0,y = —1/2, as in Figs. [2] and |4] below. We therefore proceed directly to 
the remaining cases. 

4.4 The other cases 

For all the cases with A < 1 the trajectory in velocity space goes to infinity. These 
tracks can therefore not constitute a real roller coaster with constant normal force 
in themselves. However as parts of more complicated tracks, spliced together with 
clothoids, they might be of interest. 

For Case 3, the parabolic case, we have A = 1 and 6q = it, so that (A— cos 8q) = 2. 
The variable that simplifies the integrations turns out to be, 



sm( 



1 — cos 8 
with range — oo < ui < 00. From Eqs. (l2lT) - (|24|) one finds, 



(41) 



t{u) = i(^ 2 + 3V, (42) 

s(w) = — (3w 4 + 10w 2 + 15)w, (43) 

15 

x(w) - i(^ 4 -5V, (44) 

5 

y{u) = -i(c 2 + l) 2 , (45) 

for this track, in terms of the parameter w. Fig. [2] shows the shape of this one-loop 
track. 

We now come to the two hyperbolic cases. For Case 4 we have A < 1 and (9 = 
so (A — cos#o) = (A — 1). Let, 



, / /1 + A sin(9 \ 

X = 2arctanh(y — TT ^j, (46) 

with range — oo < x < oo. Then, 

*(*> A ) = (1 ( _ X A 2)3/ 2 ( A X + Sinh x), (47) 

s(X,A) = --^L[(2A 2 + l) X +(4A + cosh X )sinh X ], (48) 
s(X.A) = ~ J A _~ A iL 2 {3Ax + [2(A 2 + 1) + Acosh x ] sinh x } , (49) 

V(X,A) = -^-|^(A + coshx) 2 . (50) 




Figure 2: Plot of the unique Case 3 trajectory with normal force equal to gravity 
(A = l). 

The shape of some of these tracks are shown in Fig. [3] One notes that in the limiting 
case of A — > this trajectory becomes, 



Kx) = - sinhx, s(x) = -(1/2) (cosh x sinhx + x)> 
x(x) = -sinhx, y(x) = -(1/2) cosh 2 x, 



(51) 
(52) 



the parabola that is the ballistic trajectory of free fall. 

One may also note that the limit A — > 1 — when taken in the x _ dependent parts 
of s(x, A) and x{x, A) make these identical: x(x, 1—) = s(x, 1 — )- This means that 
the curve becomes a straight horizontal line. 

Finally we have Case 5. Here A < 1 and 9q = it so (A — cos#o) = (A + 1). Let, 



rj = — 2 arctanh 



1 - A sin 6 
1 + A 1 - cos t 



(53) 



with range — oo < rj < oo. Then 

(A + l) 



s(t],X) 
x(r),X) 



(1 - A 2 ) 3 / 2 

(A + l) 2 
2(1 -A 2 ) 5 / 2 

(A + l) 2 
2(1 - A 2 ) 5 / 2 
(A + l) 2 



(-Xr) + swhr)) , (54) 

[(2A 2 + l)r7 " (4A - cosh n) sinh rr\ , (55) 

{3Ar/-[2(A 2 + 1) - A cosh rj] sinh r/} , (56) 

(A — cosh 77) . (57) 



2(1 -A 2 ) 2 

For A = this curve also becomes a the ballistic parabola. For positive values of A 
this trajectory has a loop and it is shown in Fig. 2] for some A- values. 

5 Conclusions 

We have presented explicit solutions for the problem of which smooth tracks produce 
a normal force of constant magnitude on a particle that slides along the track under 





Figure 3: Plot of a Case 4 trajectories. The lowest (innermost) red curve is the 
parabola of free fall corresponding to A = 0. Above that one the curves correspond 
to A = N/mg = 1/4 (green), A = 1/2 (brown), A = 3/4 (blue) and finally the 
horizontal straight line (violet) corresponds to A —¥ 1. 




Figure 4: Plot of a Case 5 trajectories. The outermost red curve is the parabola 
corresponding to N = A = 0. Inside the parabola there are curves corresponding to 
the A- values 1/4 (green), 1/2 (brown), and 3/4 (blue). The limit A ->• 1 is in Fig. [2] 



the influence of a constant gravitational field in a vertical plane. The obvious 
application is to the design of roller coasters. As we hope to have shown, the 
problem in itself has many interesting and surprising features. One of them is the 
relation to the Kepler problem of planetary motion. Students can learn various 
classical techniques of analysis by studying this problem. One notes that a direct 
numerical approach cannot reveal the fact that there are five (or four, depending 
on the point of view) qualitatively different regions of parameter values, nor the 
qualitative features of the corresponding solutions and their limiting behavior. 



A Appendix 

Explicit results for the integrals in Eqs. (f2"Tj) - (|23|) . without the constants factors 
(A — cos#o)'\ with k = 1 or 2, in front, i.e. 

d0' 



ti(9,X) = 

Si (e,x) = 
Xi(e,x) = 



9o (A-cos6") 2 ' 
„„ (A-cos6") 3 ' 

COS0W 



l eo (A-cos0') 3 ' 

are given here. For each of the three integrals there are five cases (l25j) - (l29t and 
i = 1, 2, 3, 4, 5 represents the Case number. 
One finds: 



h(e,x) 

t 2 (6,X) 



1 




(A 2 - 
1 


1) 


(A 2 - 


1) 



VA 2- 1 



... 1 (2-cos6»)sin6' 

h{0) = --- 




3 (1 — cos( 



t 4 (M) 



t 5 (e,\) = 



1 




(A 2 - 
1 


1) 


(A 2 - 


1) 



2A , / /1 + A sin( 

: arctanh 



sm 9 



Vl^A^ 

2A 
7T = A^ 



arctanh 



l-Al + cos6»y (A-cos6>) 
i sin 9 



1-A sin( 



l + Al-cos0/ (A-cos6>) 



si(6,X) 
s 2 (9,X) 

s 4 (6,X) 

s 5 (e,x) 





1 


(A 2 


-1)2 
1 


(A 2 


-1)2 



2A 2 



1 



VA2^T 

2A 2 + 1 
~\/A T— T 



: arctan 



arctan 



1 ^cos^-Gcos^ + Tlsin^ 




15 



1 




(A 2 - 
1 


1)2 


(A2- 


1)2 



(1-COS0 

2A 2 + 1 

2A 2 + 1 
"VI - A 2 " 



arctanh 



arctanh 



'1 + A 


sin 9 




1- 


-Al 


+ cos 
sin0 





'l- 


-A 





1 + A 1 - cos ( 



(4A 2 -3A cos 9- 1) sin 6» 



2(A-cos6>) 2 
(4A 2 -3Acos6»- 1) sin 9 
2(A-cos6>) 2 



(4A 2 -3A< 



1) sin( 



2(A-cos6»)^ 
(4A 2 -3Acos<9- l)sin( 



2(A-cos( 



xi(6,\) = 

x 2 (6,\) = 

xz{6) = 

Xi(6,X) = 

x 5 (6,X) = 





1 




(A 2 


1 


1)2 


(A 2 


- 


1)2 



3A 



v / A 2 ~^T 
3A 



y/W^l 



arctan 



: arctan 



A + l sin0 \ ((2A 2 + l)A-(A 2 + 2)cos6»)sin6'' 
A- 11 + costf J 2(A-cos6») 2 

', ((2A 2 + 1)A- (A 2 + 2) cos 6») sin 



A - 1 sin 9 

A + 1 1 - cos ( 



1 (cos6> 2 -3cos# + l)sin( 



(1 — cos( 





1 


(A 2 


-1)2 
1 


(A 2 


-1)2 



3A 



VT^A2 
3A 

7T = a 1 



: arctanh 



1 + A sin 9 
1 - A 1 + cos ( 



arctanh 



1 - A sin 9 
1 + A 1 - cos ( 



2(A-cos6») 



((2A 2 + 1)A - (A 2 + 2) cos 9) sin( 
2 (A- cos 9) 2 

((2A 2 + 1)A - (A 2 + 2) cos 9) sin* 
2 (A- cos 9f 



These complicated expressions simplify considerably after transformation to new 
integration variables corresponding to eccentric anomalies as discussed above. 
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